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Abstract 

We compute the top mass dependent NLO strong interaction matching con- 
ditions to the AF = 2 effective Hamiltonian in the general MSSM. We study the 
relevance of such corrections, comparing its size with that of previously known NLO 
corrections in the limit mt — > 0, in scenarios with degeneracy, alignment, and hier- 
archical squarks. We find that, while these corrections are generally small, there are 
regions in the parameter space where the contributions to the Wilson coefficients 
C\ and C4 could partially overcome the expected suppression m t /Msi/SY- 

1 Introduction 

The phenomenon of neutral meson mixing has been the subject of extensive studies within 
the field of flavor physics in the last decades, and has provided a good deal of knowledge 
concerning the flavor structure of the SM and of its extensions pQ . Experimental data on 
the mixing of K and Bd mesons constrain very strongly the existence of sources of flavor 
violation beyond Minimal Flavor Violation within the first to families, pushing the scale 
of generic flavor violating New Physics far from the TeV scale [21 El H]. 

On the other hand, the knowledge we have of flavor transitions between the second and 
the third generation of quarks is not so precise. Besides theoretical reasons to expect New 
Physics to show its face in these type of transitions, there are already several experimental 
tensions that could be the effect of physics beyond the SM (see for example [6]). An 
important focus of the LHC era is to probe B s physics to a novel level of precision, with the 
objective of clarifying these issues. In this context, the progress towards higher theoretical 
precision in B s physics is necessary. 

In the case of the Minimal Supersymmetric Standard Model (MSSM) without model 
dependent assumptions for the mecanism of supersymmetry breaking, the AF = 2 effec- 
tive Hamiltonian describing neutral meson mixing is known to NLO in strong interactions. 
Leading order (LO) strong interaction matching conditions in the MSSM have been known 
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for some time, and were computed in Refs. [7J El they arise from squark-gluino box 
diagrams. The corresponding AF = 2 next-to-leading order (NLO) corrections come from 
two loop diagrams, and have been computed in Ref. [10] within the Mass Insertion Ap- 
proximation, and in the general case in Refs. [HI [12]. The anomalous dimension matrix 
for the complete set of operators has been computed at NLO in QCD in Refs. [331 E]> 
and a resummation of large logarithms related to large squark mass ratios in scenarios 
with hierarchical squark masses can be found in Ref. [15]. A taste of the flavor con- 
straints on the SUSY parameter space in the general MSSM can be found, for example, 

in Refs. pa esi irn ubi . 

The present time is a particularly appropriate moment to reanalyze the phenomenology 
of SUSY models because of the direct searches for SUSY particles that are being performed 
at the LHC. The experimental collaborations ATLAS and CMS have already reported 
exclusion limits based on ~ 1/fb of data in searches for jets + I$t with zero [191 12H] or 
more charged leptons, and larger datasets up to ~ 15/fb are expected in the future. 

These limits, however, are model dependent, and the lower bounds on the SUSY masses 
vary from model to model. In the constrained MSSM (CMSSM), squarks and gluinos of 
equal mass are excluded below ~ 950 GeV [19] , but the bound on the gluino mass would 
be slightly higher (~ 1.1 TeV) for squarks around 500 GeV [20]. A more detailed analysis 
including LHCb data, as well as electroweak precision and I?-physics observables obtains 
best fit values of (m 1/2 ,m ) ~ (780,450) for the CMSSM and (m 1/2 ,m ) ~ (730, 150) for 
the NUHM1 [21]. In the case of the phenomenological MSSM (pMSSM), detailed analyses 
find that 85% of the parameter space is excluded for gluino masses below ~ 520 GeV [22] . 
while rrig g > 1.1 TeV if m g ~ rriq and m g > 700 GeV if rriq ^> m g [23]. Thus, the SUSY 
scale can still be low enough for flavor physics to impose substantial constraints on flavor- 
violating soft supersymmetry-breaking couplings in general versions of the MSSM, and 
these constraints should certainly be reevaluated. At the same time this reevaluation 
should involve, as argued above, a novel level of theoretical precision. 

In Ref. [11] the AF = 2 NLO strong interaction matching conditions were computed, 
putting all the quark masses to zero. This approximation is entirely justified since correc- 
tions to this limit are suppressed by a ratio m q / Msusy *C 1. In the case of the top quark 
this suppression still exists, although it might not be so effective if the SUSY scale is 
below the TeV, with m t /MsusY ~ 1/6. In such a case, it may be that the loop functions 
have numerical enhancements that compensate such suppression: a factor of ~ 10 is not 
unreasonable, at least in some regions of the parameter space. The only way to know if 
this happens is to compute explicitly these corrections. 

The purpose of this paper is to clarify this issue, by making the two loop calculation 
and comparing the size of the corrections with the m t — > limit. After reviewing briefly 
the structure of the AF = 2 effective Hamiltonian in Section [2, we will describe the 
computation of the m r dependent corrections in Section [31 specifying the structure of the 
Wilson coefficients. The explicit expressions for the m t -dependent Wilson coefficients in 
the general case are presented in Appendix rAl and in the Mass Insertion Approximation 
in Appendix [Bj In Section 0] we analyze the size of the corrections in three different SUSY 
scenarios: degeneracy, alignment, and a model with hierarchical sfermions. Conclusions 
and a summary of the results are presented in Section [5j 
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2 Effective Hamiltonian for AF = 2 processes 



The most general effective Hamiltonian for AF = 2 processes up to operators of dimension 
six can be written as 53 

i=l i=l 

where Cj are the Wilson coefficients and Oi are the dimension six AF = 2 operators. In 
four dimensions there are eight independent operators of this type. Here we choose the 
following basis: 

01 = s a -f^P L b a s^P L b p 61 = s a ^^P R b a sp^P R bfi 

2 = s a P L b a s p P L bp & 2 = s a P R b a s p P R bp 

(2) 

3 = SaPibp Sf}P L b a e> 3 = s a P R bp s p P R b a 

O a = s a P L b a s p P R bp 5 = s a P L b p s p P R b a 

where Pl,r = (lT7s)/2 are the usual chiral projectors. To simplify the notation through- 
out the paper we focus on the case of B s — B s mixing, while other cases of interest can 
be recovered by obvious substitutions of quark fields. 

The mixing amplitude is obtained by taking the matrix element of the effective Hamil- 
tonian between B and B states: 

A = = E C ^ U ( m ^ A )ij (Oi(m b )) (3) 

i ij 

where A is the matching scale at which the Wilson coefficients are computed. The matrix 
U(rrib, A) is the evolution matrix, that accounts for the renormalization group running of 
the coefficients down to the hadronic scale. The evolution matrix at NLO can be found 
for example, in Ref. [23]. The matrix elements of the operators can be found in Ref. [2]. 
These matrix elements are calculated in the RI-MOM scheme, which means that the 
Wilson coefficients used in Eq. must be computed in the this scheme. 

The NLO Wilson coefficients presented here and in Ref. [11] are given in the NDR 
scheme. The prescription to translate between NDR, DRED and RI-MOM schemes in- 
volves rotating the set of coefficients with appropriate 0(a s ) matrices, which are given 
for example in Ref. [TO] . 

3 Top mass dependent corrections 

At two loops, the top quark appears for the first time inside the Feynman diagrams that 
contribute to the AF = 2 effective Hamiltonian. It appears in the one loop correction to 
the gluino propagator, in box diagrams such as those depicted in Fig. [TJ 
A diagram of this kind is proportional to the following two loop integral: 

^ f ,4 ,4 [^1^ + m § )t/ fc *(A^ + m t )y fc %+m § )C/*]®[^ s (^ + m § )^ 6 ] ... 
V ~ d qi dq 2 7-2 ovw , ( 4 ) 



(qf — m^) 3 (ql — mf){q\ — m 2 -)^ — m 2 k )(Aq 2 — m 
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Figure 1: A representative sample of the two loop diagrams that contribute to the 0(oP s ) 
top mass dependent corrections to the mixing of B s mesons. 



where Aq = q x — q 2 , U tq = T l ^P L — T l R P R , and V tq = T l ^*P R — T l ^*P L . The squark rotation 
matrices T are defined as being the rotations that relate squark fields in the super-CKM 
basis (qf L , q[ R ) to the mass eigenstates qf 

a i,L — 1 D L U J i u i,R ~ 1 D R U 3 

u{l = T u* L Uj , u[ R = Tj}*^ (5) 

The masses and momenta can be normalized to the gluino mass: defining Xj = m^/m? 
and x t = m 2 /m~, we have 

v ~ — ! d%d% [V1S{ ^ + ^HHHM + V^)v k \ji + i)u ib ] ® [v^i + i)u jb ] (6) 

mj J 12 (qj-l) 3 (qf- Xi)(qf- Xj)(qi- x k )(Aq 2 - x t ) 



The leading corrections for xt — > where computed in Ref. [IT]. Expanding up to 

C(m^/m?) one gets 

^ ^ \/%~t f l4 ■■■ U kt V kt ■■■ x t f j4 • ■ ■ U kt Ai V kt ■■■ 

» ~ + ^ / a„ + ^ y am ( ...^ 4 + • • ■ (7) 

where we have kept the pieces from the squark-top loop. 

We see that the mt/m g contribution is proportional to r^T^* and r^T^ , while the 
m t/ m l contribution is proportional to p^T^ and r^Tj^. This means that the m t jm g 
correction is suppressed by a LR mass insertion, as expected for the required chiral flip. 

The corresponding contributions to the Wilson coefficients in Eq. ([1]) from all these 
diagrams can be written as 



3 r 

/nr(t) _ a s ijkf-pijk , -pijk \ , m t i ijk-pijk , i ijk-pijk 

~ 2 / j a £ V ABC DLL ' ABCDRR/ rn _ ABCDLR ' £ AB 



12%m? 

9 ijk 



m-g 



ABC DLL) 



J _ r Jh_fijk/ v ijk j_-pijk \ , 

" t " rrl 2 h v 1 ABCDLL 1 ABCDRR) 



where ^ % abcdef — ^a*^b^c*^d^e*^f (capital letters mean L or R), and the sum runs 
over 6 down-type squarks in the case of i,j and 6 up-type squarks for k. 
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In Eq. (jSJ), the functions af correspond to the m t = contribution, and are given 
explicitly in Ref. [TT]. The functions f^ k , gf k and h^ k have been calculated here, and 
are given explicitly in Appendix [A] The details of this calculation are the same as those 
described in Ref. [TT], and we refer the interested reader to that reference. 



4 Size of the corrections 

In this section we will study the size of the m r dependent corrections. Naively, one expects 
these corrections to be at the few percent level with respect to the m t = correction: 
the first term is of order mt/rrig < 20%, but contains a further suppression from the LR 
mass insertion; the second term is of order m^/m- < few%. Constraints on the LR mass 
insertion come from vacuum stability considerations, which require that [25] : 

^<2f, (9) 

but the exact bounds depend on the particular spectrum of the model. Therefore, the 
two contributions have roughly the same suppression. The question is whether the loop 
functions can account for a compensating enhancement, at least over some regions of the 
parameter space. 

In order to address this issue, we will explore the SUSY parameter space in three 
different scenarios: degeneracy, alignment and hierarchy. 



4.1 Degeneracy 

We consider the Wilson coefficients in the Mass Insertion Approximation, where the 
squark squared mass matrix in the super-CKM basis has degenerate diagonal entries much 
bigger than the off diagonal terms (mass insertions). The Wilson coefficients depend on 
a reduced set of parameters: the average squark mass M s , the gluino mass — through 
the parameter x = m~/M^, and the mass insertions 5^ L , S R b R , S rl , S^ r and 8 l l R . The 
Wilson coefficients in the MIA can be obtained from the general expressions presented in 
Appendix [A] The general form of these coefficients is: 



C 



MIA 



a 



3 



12 3 ttM? 



2 

(H{x) + — gi (x) Re{Sf R ) + ^| fi(x) + ■ ■ 
rrig m~ 



Md , ( 10 ) 



where A, B,C,D = L, R. The functions aj(x) correspond to the m t = contribution, and 
can be found in Ref. [T0]|3 The functions fi(x) and gi(x) can be found in Appendix IB"1 
where the complete expressions for the Wilson coefficients in the MIA are given. 

In order to estimate the size of the new corrections, we consider the following ratios, 
Rf(x) = gi{x)/ai{x) and R{(x) = fi(x)/a,i(x), and determine if there is a region in the 
variable x where any of these ratios are large. 



1 In this reference the coefficients are given in the DRED scheme. One should convert them to the 
NDR scheme to be consistent with the conventions in the present paper. 
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Figure 2: Relative size of the m^-dependent contributions in the Mass Insertion Approximation. 
The ratios R\ are shown as solid lines and R{ as dashed lines. In the case of R®, the upper 
solid curve shows the relative size of the contribution multiplying S^ L and the lower solid line 
the one multiplying 5^. In the case of coefficients C4 and C5, the two possibilities (LL + RR 
or RL + LR insertions) are displayed. 



In Fig. |2]we show the Ratios Rf and R{ as a function of x for i = 1, .., 5. We see that 
the only contribution that can partially overcome the suppression is the contribution to 
C\. In the region x ~ [0.5 — 1], this contribution could be up to ~ 20% of the whole NLO 
contribution, if m- g ~ 1 TeV and ~ 0.2. The same is obviously true for the coefficient 
C\. The rest of the contributions can be safely neglected over the whole parameter space 
in the degenerate scenario. 
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4.2 Alignment 

A second scenario to consider is the "aligned" scenario, in which the squark mass eigenba- 
sis is approximately aligned with the super-CKM basis. In this case, the diagonal squark 
masses can be different, with wide mass splittings. However, since the mixing angles 
are small (the rotation matrices are close to the identity), the off-diagonal entries (mass 
insertions) of the squark mass matrix in the super-CKM basis are small. An expansion 
in these mass insertions can be called the Non-degenerate MIA, or NDMIA. 

From the full results of the Wilson coefficients presented in Appendix |A] and in 
Ref. [TT], the expressions of the WC's in the NDMIA can be obtained (see for exam- 
ple Section 6 of Ref. [TT]). These coefficients have the same structure as the MIA WC's, 
but they depend on a set of different squark masses instead of on a single variable x: 



£yNDMIA _ a "s 



12 3 7rm? 



a,(x) + ^ 9i (y) Re(6* R ) + 4 /i(y) + ■ ■ 
m g m~ 



S AB S CD » i 11 ) 



where the m t = part depends on all 12 squark masses (as they all appear in the 
squark-quark loop): x = (xq L ,Xg R ), and the m t ^ part depends on at most 6 masses: 
y = ( x s A , ^i B , x Sc , xi D , xi L ). Also, the barred mass insertions 5 l ^ Y are the mass insertions 
normalized to the gluino mass (instead of to an average squark mass, which is now ill 
defined) . 

The large number of free parameters in this scenario does not allow for a consistent 
scan of the full parameter space. I any case, we are interested in finding "tipical" regions 
in which the new contributions are not negligible, so a full scan is not necessary. Since 
we are focusing on B s mixing, the most relevant squark masses are m s and mj, so the 
strategy is the following. We scan randomly over squark masses (up to several times the 
gluino mass) for all squark masses except m s and m^, and we plot in the m~ a -m~ h plane the 
ratios R\ = gi/ai and R{ = jija^ 

In the left plot of Fig. [3] we show the ratio Rf as a function of and for a 
random set of squark masses: xq l = 1.2, xn R = 0.6, xg L = 0.7, x^ R = 1.8, xs R = 0.2, 
Xc L = 1.9, Xc R = 0.8, x~ bR = 0.4, x~ tL = 0.3 xi R = 0.5. The dashed line bounds the region 
where Rf > 2 and the thick solid line indicates the region where R\ > 3. In this region the 
size of the new contribution can be as high as ~ 10-15% of the whole NLO contribution to 
the Wilson coefficient C\. The same is true again for C\. The ratio R{, on the contrary, 
is in general small and can be neglected. 

The contributions to the other Wilson coefficients are all negligible, or at least we 
have not found any generic region where the ratios Ri for % > 2 are bigger than 1. As an 
example we show, in the right panel of Fig. [31 the ratio i?f in the case of LR, RL mixing. 
The squark masses used are the same as in the previous case. We see that this ratio is 
tipically smaller than 1, and around ~ 0.5 for ms,m~ b ~ 0(rrig). Very similar results are 
obtained for the rest of the coefficients. 
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Figure 3: Relative size of the mi-dependent contributions in the Non-degenerate Mass Insertion 
Approximation. Left: Rf as a function of Xg L and x-^ . The dashed line bounds the region Rf > 2 
and the thick solid line bounds the region Rf > 3. Right: i?| in the case of LR, RL mixing. In 
both cases we take the rest of the squark masses ~ O(m^), as indicated in the text. 

4.3 Hierarchical squark masses 

A third scenario to consider is the case of hierarchical squark masses. In this scenario, 
the first two generations of squarks are assumed to be much heavier than the rest of the 
supersymmetric spectrum, lying near the TeV scale. Besides satisfying naturalness criteria 
[26] . this scenario avoids flavor bounds from processes involving the first two families. 
This is due to the fact that contributions mediated by first two generation squarks are 
suppressed by a small ratio of masses, and the processes must proceed through the third 
generation, at the cost of containing an extra mass insertion. 

A convenient parameterization of the squark rotation matrices in this scenario is the 
following [TB~] : 





5l 3 £cos6 



1 

5fi cos e 


°LL 

°LL 
COS 9 





<5££ sin 9e^ 





<5f^ sin 9e i(t> 





sin 9e 1 ^ 









-5f* sin Be-** 





— sin 9e~^ 


1 


5 l ^cos9 



1 

^rr cos 9 


°RR 

_A23 
°RR 

COS 9 



(12) 

where it should be understood that there are two rotation matrices (IV and r^), each 

8 




4 6 8 10 12 14 1.5 2.0 2.5 3.0 3.5 2 4 6 8 10 12 14 

x e x £ xi 



Figure 4: Relative size of the mj-dependent contributions to C4 in the Hierarchical squark mass 
scenario: i?| as a function of xi and Xf R . We show an overview plot of a large region (left), 
and two particular regions where the ratio is large and positive (center) and large and negative 
(right). In these plots we have taken Xh = 1000, 9 = 45° and eft = 0. 

with its own set of parameters. The 8's are dimensionless mass insertions, and 9jj,d and 
4>u,d are mixing angles and phases that determine the diagonalization of the light-squark 
sector. The parameters involved in this scenario are: a common squark mass for the heavy 
sector fhh ~ 10 — 100 TeV, a common squark mass for the light sector rri£ ~ 1 TeV, two 
mixing angles 9u,d and two phases 4>u,d describing the mixing in the light sector, and the 
mass insertions Sf LRR , S s L b LRR , 8f LRR and 5f LRR . 

In this particular scenario, with this parameterization, and with no mass splitting 
between third generation w-type squarks, the corrections proportional to m t /mg are zero: 

#(zfc) r L* r K = g{x iL ) sin9 cos9e^ - g{x iR ) sin9 cos9e^ = . (13) 

k 

The loop functions fi related to the corrections proportional to mf/m- are, as happens in 
the other scenarios discussed previously, small with respect to their m t — > counterparts. 
We have checked that the ratios R{ in this case are below 10% in all cases. 

A mass splitting between ti and t R can change the situation. We choose the "right 
handed" stop to be the one whose mass deviates from ihg. In this case the ratios Rf are 
nonzero for Ci, C4 and C5. We find that these ratios are not more than a few percent for 
C\ and C5, but they can be moderately large in the case of C4. In Fig. H]we show the 
values of R\ as a function of X£ and taking for illustration for the other parameters 
the following values: Xh = 1000, 9 = 45° and = 0. While this ratio is zero for fhp = m^, 
it increases with the mass splitting, ranging from ~ —4 to ~ 10 in a relatively wide region. 

5 Conclusions 

In this paper we have computed the top mass dependent two-loop squark-gluino correc- 
tions to the AF = 2 matching conditions. The expanded results for the Wilson coefficients 
up to order m^/m^ are given explicitly in Appendix El (for the general case), and in Ap- 
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pendix[B]in the Mass Insertion Approximation. This complements the results of Ref. [TT] . 
that were obtained in the limit m t /m g — \ 0. 

We have studied the relative size of this corrections compared to the m t — > results. 
For that matter we have considered three different SUSY scenarios. In the case of degen- 
eracy, we find that the corrections are small and can be neglected in general, but there 
is a region for m 2 g /M 2 ~ 0.5 — 0.9 where the contribution to the coefficient C\ could be 
notable. 

In the case of alignment, the situation is very similar. All contributions are negligible, 
but there is again a region in which the contribution to C\ could partially overcome the 
suppression, corresponding to values of the squark masses Xs L ,x~ bL ~ 0.5 — 1.5. This is 
true for a chosen set of values for the rest of the squark masses, which are specified in the 
text. 

In the case of the hierarchical squark mass scenario, the situation is quite different. The 
correction linear in m t /m g is strictly zero for degenerate stop quarks, given the chosen 
parameterization of the squark rotation matrices. When we allow for a mass splitting 
between left and right handed stops, these contributions are nonzero. We find that in this 
case they are negligible for all the coefficients except for C4. The contribution to C4 can 
become relevant in some regions of the parameter space specified in Fig. HI 
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A Wilson Coefficients 

In this section we give the explicit expressions for the Wilson coefficients. We omit the 
results for 6*1,2,3, which are obtained from 6*1,2,3 by exchanging L -H- R in the squark 
rotation matrices. Here we give only the m t dependent part of the Wilson coefficients: 

C™ - * = Ci — Ci\ mt= o (14) 

The functions depend on the squark masses through the squared mass ratios Xi , defined 
aS Xi = TtX,: /m~. A sum over k is understood, running over all right and left-handed 
squarks of u or d type depending on the case. We also define the product of squark 
rotation matrices: 

-pijk = -pis*-pib -pjs*-pjb -pkt*-pkt /ir\ 

1 ABCDEF — 1 A B C D E 1 F 

Expanding the Wilson coefficients in powers of m t /mg we have, up to order (m t /m g ) 2 : 



C\ 



a" 



l2Hm 2 g 



m t 



rrih 



gf (r 



LLLLLR 



_|_ pijfe 



m 



LLLLRLj 



m 



ZT2 ■/! V LLLLLL ' 1 



LLLLRR) 



l2Hm 2 g 



rrih 



ijk -pijk , 7 ijk -pijk \ 

h 1 RLRLLR i" n 2 1 RLRLRL J 
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m t rijk fpijk 

2^2 \ RLRLLL 

nig 



+ ^ rlrlrr) + 



(16) 
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-yijk 
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^ijk > 
LRRLRR, 
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_|_ 

The loop functions can be written as 
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(17) 



All these can be expressed in terms of 5 master functions, that we have chosen to be Fi, 
F 2 , Gi, G 2 and A, such that: 
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The master functions are given by: 
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x fc log x fc 



(20) 



+ 



19x? - 57x fc + 6* 



2 i 2 

x fc log X k - 



( x . _ l)3( Xfc _ 1)3 

2x 4 (llx 4 + 19) 



2x fc ((19x fc - 15)4 - ( 4 9^fc - 26)xi + 19) 



logx fc 



fa - l) 3 



log Xi - 2 



( x . - l)2( Xfe _ 1)2 

{llx k + 4)4 - {llxj - 15x fc + 34)xj - x k (19x k - 34) 
(x,-l) 2 (x fc -l) 2 



F 2 ife = 816 



2xiX fc (x 2 - 3xi + 3) 



fa - l) 3 



Li 2 (l - x k ) 



- 2x " Li 2 (i - a; 



X3 o 2 i o , _ , > • 3 i o ,->-• 2 o ™» . 

fa - l) 3 fa - l) 3 



x fc log 2 Xfc 



x fc (xj - 2)(2x k Xj - Xj-l) 
fa - l) 2 (x fc - l) 2 



logx fc 



2xi xf + fa - 3fa - 2x^ + 3x fc 

log - - 



fa - l) 3 



(x,-l) 2 (xfc-l) 2 



(21) 



Gf = 24 



2x,(-2x fc x 2 + (6x fc + 13fa - 19x fc + 2) 

Li 2 (l - x k ) 



fa - If 
2(13x i + 2)(x l -Xfc) _ . 



fa - 1) £ 



_ 2x 3 k ~ 6x 2 k + 4 5zfc + 4 

121 xk) fa - i) 3 fa - 1) 3 



x k \og x k 



3(13xfc + 2)x k Xi + 13x^ — 156x 2 + 21xfc — 13 2 2 Xi(13xj + 2) 



fa - i) 3 fa - ly 



Xi log x fc - 



(xi - iy 



log 2 Xj 



13xi - 41x2 - 72x? - 37x fe + 2 , 2 2x^13^ + 2), 

Xi log Xfc / -T^— log Xi log Xfc 



(Xi - l)3(x fc - 1)3 & (Xi-1) 3 

(2x| - 6x 2 k + 110x fe - 5)xi - 2xf - 9x| - 80x fe - 10 



(x,-l)(x fc -l) 3 
(219x fc - 17)xj - {219x 2 k + 185 fa + x fc (17x fc + 185) 



. , Xi(89xi + 12) 
x fc log Xfc H — - — log Xi 



fa - l) 3 



2(x 4 - l) 2 (xfc - l) 2 



(22) 



Gi k = 816 



2xi(xi Xfc) . , . 2xj(xi Xfc) . , . 
- 1)3 2 " ^ + (X, - 1). " £) 



Q ry»2 ™,3 | ( 3 1 Om2 I Q ,-y> __ T A ,y»2 / ry*4 Q / >-» 3 ft ^y.2 ,-y» , \ ,->-» . Q 'v>2 

OJsfeU;^ ~r \ ' k ) i \ k k k k 1 ox^, 

(x,-l)3(Xfc-l)3 

2x? (7xfc - l)xi - x 2 - 5xfc 

log Xi log Xfc — — Xfc log Xfc 



x- 



{xi - ly 



(Xi-l)(Xfc-l) 3 



(Xi - l) 3 



log 2 Xfc 



log 2 Xi 
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6a; 2 . {7x k - l)x 2 - (7x 2 k + 5)xi + x k (x k + 5) 
log Xi H 



(x,-i) 3 "°- ' (xi - i) 2 (z fc - 1; 

Here, the function Li 2 (x) denotes the dilogarithm defined in the usual way: 

log(l-t) 



Li 2 (x) 



dt 



(23) 



(24) 



B Wilson Coefficients in the MIA 

In this Section we present the expressions for the Wilson coefficients in the Mass Insertion 
Approximation. Again, we present only the m t dependent part of the Wilson coefficients: 



(J mt 



i \mt=0 



(25) 



The coefficients Ci.2,3 are obtained from by exchanging Lofiin the squark rotation 
matrices, and C™* = —3/17C™*. The variable M 2 denotes the squared squark mass, 
meaning any of the common diagonal entries in the squark mass matrix in the super 
CKM basis (the diagonal masses are corrected by the mass insertions). The coefficient 
functions depend on the variable x = rn^/M 2 . The rest of the Wilson coefficients are 
given by: 



ntt 



a 3 m t 



12 3 ttM 2 m~ g 3(x - l) 6 



a 



nit 



11360a; + 2056a; 2 - 14496x 3 + 1144a; 4 - 64a; 5 

- 3744x logx - 15696x 2 logz - 14976xLi 2 (l - x) - 28224x 2 Li 2 (l - xfjRe(6^ R ) 

+ — ( - 10208a; - 44032x 2 + 51576a; 3 + 2816a; 4 - 152a; 5 - 2592a; 2 log x 
m- g \ 

- 21168a; 3 log a; + 6336xLi 2 (l - x) + 36864x 2 Li 2 (l - x 
1 



(26) 



av 



m t 



(^15232a; - 2856a; 2 - 13056a; 3 + 680a; 4 - 4896a; log x 



' \2HM 2 m g 3(x- l) 6 

- 17136a; 2 log x - 19584a;Li 2 (l - x) - 29376x 2 Li 2 (l - x])5% R + ( - 2448; 
+ 36040a; 2 - 38352a; 3 + 5304a; 4 - 544a; 5 - 22032a; 2 log a; - 48960x 2 Li 2 (l - xj)6%£ 



+ —(- 143888a; 2 + 135456a; 3 + 8976a; 4 - 544x 5 + 4896a; 2 log x - 58752a; 3 log x 



CI 



+ 97920a; 2 Li 2 (l -x 



12 3 ttM 2 ^g 3{x - 1] 



«) 2 (27) 
20832a; + 84000a; 2 - 117792a; 3 + 14304a; 4 - 1344a; 5 
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- 8640a; log a; - 84672a; 2 log a; - 34560xLi 2 (l - x) - 172800x 2 Li 2 (l - x)^Re(5 
+ — flll36z - 325920a; 2 + 295776a; 3 + 20256a; 4 - 1248a; 5 + 15552a; 2 log a; 



tt \ 

LR) 



- 129600a; 3 log x - 6912xLi 2 (l - x) + 214272a; 2 Li 2 (l - x] 



+ 19712a; + 3696a; 2 + 16896a: 3 - 880a; 4 + 6336a; log x + 22176a; 2 log x 

+ 25344xLi 2 (l - x) + 38016x 2 Li 2 (l - x)^Re(^ R ) + ^(^20416x + 54208a; 2 
- 71280a; 3 - 3520a; 4 + 176x 5 + 6336a; 2 logx + 28512a; 3 log a; - 12672xLi 2 (l - x) 



a 



50688x 2 Li 2 (l - x) 
a 3 m t 1 



sb \(xsb > 
RLj 



«r)(5 



(2? 



19424a; + 544a; 2 - 21408a; 3 + 1504a; 4 - 64x 5 



" 12 3 ttM s 2 m § 3{x - l) 6 

- 6336a; log a; - 24768a; 2 log a; - 25344xLi 2 (l - x) - 43776x 2 Li 2 (l - x)\Re{S^ R ) 
+ —( - 18560a; - 66208a; 2 + 80736a; 3 + 4256x 4 - 224x 5 - 5184x 2 log x 



32832a; 3 log a; + 11520xLi 2 (l - x) + 57600x 2 Li 2 (l - x) 



{^llWrr) 



- 26880x + 5040a; 2 + 23040a; 3 - 1200a; 4 + 8640a; log a; + 30240a; 2 log x 



+ 34560xLi 2 (l - x) + 51840x 2 Li 2 (l - x)]Re(5f R ) + — (27840a; + 73920a; 2 

/ nig V 

- 97200a; 3 - 4800a; 4 + 240a; 5 + 8640a; 2 log x + 38880a; 3 log x - 17280xLi 2 (l - x) 



69120x 2 Li 2 (l -x 



{^lrWrl) 



(29) 
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